Spatial correlations in finite samples revealed by Coulomb explosion 
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We demonstrate that fast removal of many electrons uncovers initial correlations of atoms in a 
finite sample through a pronounced peak in the kinetic-energy spectrum of the exploding ions. This 
maximum is the result of an intricate interplay between the composition of the system from discrete 
particles and its boundary. The formation of the peak can be described analytically, allowing for 
correlations beyond a mean-field reference model. It can be experimentally detected with short and 
intense light pulses from 4th- generation light sources. 
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Correlation effects beyond the mean-field description 
are important in understanding the detailed structure 
of matter. A major scheme has been the "exchange- 
correlation hole" indispensable in electronic-structure 
theory [1 . While this research focuses on bound, mostly 
even on the ground state, very intense light pulses now 
available for XUV to X-ray energies [2 , create situations 
to probe correlation in the continuum. One phenomenon 
in this context is massively parallel ionization [3 of a 
cluster or large molecule: Many electrons are ionized al- 
most simultaneously to an energy high enough such that 
the attractive ionic background does not play a role on 
the short time scale of energy exchange and correlation 
of electrons in the continuum. 

Complement arily, one may ask how the many ions be- 
have whose dynamics develops on a longer time scale 
when most of the electrons have left. Do the ions also 
develop a correlated motion in the continuum — and if 
so, what is its characteristics in terms of observables, 
e.g., the energy spectrum of the ions? The simplest de- 
scription well-known from plasma physics approximates 
these ions as a continuous charge density. The corre- 
sponding spectrum provides a comparison for highlight- 
ing the diflFerences to a correlated continuum calculation. 
They occur for two reasons: Firstly, the granularity of 
our system, i.e., the fact that we deal with real particles 
which typically have a minimal separation 5 and therefore 
give rise to a classical "correlation hole". Moreover, the 
ionic system is finite (for simplicity we assume a spheri- 
cal shape) with an edge. The softness of this edge, i.e., 
the distance a over which the charge density decreases to 
zero, is the second spatial scale which influences the ionic 
correlations in the continuum. In fact, we will show that 
the ion dynamics exhibits a crossover, characterized by 
the ratio a/5, from a dominant influence of the correla- 
tion hole if a <C (hard edge) to the dominant influence 
of the edge for 5 (soft edge). 

Surprisingly, this crossover can be probed by simply 
varying the pulse length of the ionizing light pulse in the 



correct parameter regime: A long pulse (for reasonable 
target systems this would be of the order of a few hun- 
dred femtoseconds) will slowly ionize the cluster and give 
the ions created near the surface time to move outward. 
Thereby, a soft edge in terms of decreasing density of ions 
is created. On the other hand, a short violent pulse (of 
a few femtoseconds only) will start the motion of all ions 
created at once leading to correlation-hole dominated dy- 
namics. 

We set the stage by deflning the mean-fleld reference 
system, a homogeneously charged sphere of charge Q and 
mass M. It undergoes Coulomb explosion, doubling its 
initial radius R in time 
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Making use of Gauss' law and the fact that shells of 
charge in a homogeneously charged sphere do not over- 
take each other, one can directly map the initial potential 
energy due to the enclosed charge at a radius r to the flnal 
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FIG. 1: (Color online) Energy spectra as obtained from an 
exploding LJ cluster with 1000 particles for various pulse du- 
rations T given in units of the expansion time t^, cf. Eq. ([T]). 
Energy is given in terms of the highest energy for ho- 
mogeneously charged sphere, cf. Eq. (3b). The dotted line 
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marks the surface energy of a sphere, which is charged homo- 
geneously according to a Gaussian pulse (see text). 
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FIG. 2: (Color online) The maximal kinetic energy of the ions 
in Fig.[l]for different pulse lengths T/t^ (a). The two arrows 
indicate the pulse lengths T — t^/20 and T — 20t^ for which 
the evolution of the spectra dP/dE are shown in (b) and (c) 
respectively as a function of energy and time ^ G [0, 1]. The 
time variable ^ = E{t)/E{oo) describes how far the Coulomb 
explosion has evolved in terms of the kinetic energy E{t) of 
an expanding homogeneously charged sphere. 



kinetic energy E per unit charge E{r) = Qr^ /R^. Evalu- 
ating the standard expression for the energy distribution 
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In terms of these scales for energy E^ and time t^, we 
show in Fig.[T] kinetic-energy spectra of Lennard- Jones 
(LJ) clusters [11] with 1000 particles for various effec- 
tive pulse duration T. Hereby, effective pulse duration 
refers to the time in which the particles get charged. 
We select randomly a charging time for each particle in 
such a way that the ensemble averaged charge grows as 
{Qt) = Q[l+erf(t/T)]/2 which, e.g., corresponds to ion- 
ization of atoms in a cluster by single-photon ionization, 
as it would occur for a sample exposed to intense free- 
electron laser radiation in a Gaussian pulse [7 . 

One clearly sees in Fig.[l] a cross-over behavior with 
two "hot spots" at short and long pulses (cf. also Fig.[2|i 
which shows the height of the ridge as a function of pulse 
length T), always close to the maximal energy possible 
at a given pulse length T. The long-pulse regime is 
compatible with a mean-field dominated dynamics and 
has been described previously: This mean-field peak in 
the kinetic-energy spectrum was observed for spherical 
objects with a washed-out edge, i.e., a non-uniform ra- 
dial density distribution. It allows faster inner regions of 
exploding charge density to overtake slower outer ones, 
leading to radial caustics, as revealed by means of a ki- 
netic model Under suitable conditions this caustic 
can even lead to a shock shell [5]. Pre-exploding the 
sample [9 softens the edge of the system to a gradually 
decreasing density, which is indeed close to the situation 



for long pulses shown in Fig.[T] Charging at low densities 
as effected by a long laser pulse reduces the energy de- 
position in the sample and leads to a global decrease of 
kinetic energies for long pulses, also visible in the spec- 
trum for a homogeneously charged sphere (dashed line in 
Fig.[l]). 

What has not been seen before is the sharp maximum 
for short pulses whose origin is fundamentally different 
from the mean- field peak at long pulses. We will refer 
to this maximum as the granularity peak since its origin 
is the granularity of the system in connection with the 
fact that it is finite and has an edge. Indeed both fea- 
tures, the mean- field and the granularity peak are sur- 
face phenomena. This fact can be easily read off from 
Fig.[l] since for the respective pulse length the peaks are 
located near the maximal energy contributed by surface 
ions. This is further confirmed by comparing the mean- 
field spectra according to Eq. Q and the numerical ones 
from individual ions [12 in Fig.|3^-c. Formally, these 
spectra correspond to cuts at T = in Fig.[l] but for dif- 
ferent cluster sizes, namely 100, 1000 and 10000 atoms, 
respectively. One sees, that the spectrum deviates less 
from the mean-field prediction for larger cluster size and 
correspondingly smaller surface-to-volume ratio, never- 
theless the granularity peak remains clearly visible. 
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FIG. 3: (Color online) Spectra from numerical propagation 
of systems with N — 100, 1000 and 10000 ions, respectively, 
are shown by orange/gray- shaded areas (a-c). The mean- 
field model as given by Eq. (|3| is shown with dashed lines. 
Considering a correlation hole in the mean- field description by 
using Eq. ^ in Eq. ([2| yields the spectrum shown as solid line 
(b). Additionally the pair-correlation function g{r) according 
to Eq. Q is shown for the cluster with A^=1000 individual 
ions (d). 
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That mean-field and granularity peak are of different 
origin can be seen, e.g., from the formation of the max- 
ima in time: The mean- field peak builds up gradually 
in time and appears towards the end of the explosion 
(Fig.|2]3) since the inner density needs time to overtake 
the outer one. The granularity peak in case of a short 
pulse, on the other hand, appears right from the begin- 
ning (Fig.[2]3) suggesting that it can be understood from 
the forces acting initially on the individual ions. This is 
indeed the case as we will see. 

Due to the repulsive interaction potential of the atoms 
in the cluster, nearby particles are very unlikely in the 
ground state. Hence, a correlation hole around an atom 
or ion exists. This means if there is an atom at po- 
sition r the probability to find another atom at posi- 
tion within the radius S of the correlation hole, i.e., 
for |r— r'l < 5, vanishes leading to the characteristic pair- 
correlation function as shown in Fig.|3]i. It has been nu- 
merically obtained by an ensemble average over 
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with Vij the distance between the particles i and j and 
g(r) = ^(r/i?-2)^(r/i?+4) the distribution of distances 
r in a sphere with radius R \W . The normalization with 
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FIG. 4: (Color online) Sketch of the integration assuming a 
correlation hole with radius 5 around the test particle at dis- 
tance r. We show two situations, where the correlation hole is 
inside the bulk (ri) and at the surface (r2), respectively. The 
test particle force is obtained by integration over all shells (in- 
dicated by green/gray lines) with some of them being "open" 
(thick green/gray line, see also inset) due to the correlation 
hole. The force on a test particle as a function of the distance 
r (a) due to an "open" charged spherical shell with radius r 
according to Eq. (|6| and (b) integrated over all shells of the 
charged sphere of radius R according to Eq. ([7|. 



g{r) is necessary in order to remove the trivial depen- 
dence on r due to the finite size of the system. 

We can assess how the correlation hole affects the force 
acting on an ion at position r analytically. To this end 
we determine the Coulomb repulsion of the ion from an 
infinitesimally thin spherical shell of radius r' and charge 
q by integrating over all angles 
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with = r'(sin6>cos^, sin6>sin^, cos6>) a vector on the 
shell. The radial component of this force can be deter- 
mined without loss of generality by choosing r = rz along 
the 2:-axis. Integratig over (j) and using r = cos it reads 
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Hereby, the integration was restricted to the upper limit 
1—St<1 in order to account for a correlation hole, cf. 
the sketch in Fig. 4] and its inset, which shows St explic- 
itly. Equation (p) corresponds to Sr = 0. In this case, 
performing the integration over r yields — as expected — 
Gauss' law, i.e. /r',o = for r < and /r',o = q/r'^ for 
r >r' . As can be seen in Fig.|4^, for finite 5t the force is 
fr',5T > inside the shell with radius r' and fr'^Sr < q/'^r'^ 
outside this shell. Yet, this modification of the forces does 
not play a role as long as the test particle's correlation 
hole is in the bulk {r < R — see test particle at ri in 
Fig.[4| since reduced repulsion from inner charged shells 
with < r is fully compensated by a finite repulsion 
from outer shells with > r in accordance with Gauss' 
law for spherical charge distributions and a test particle 
with its correlation sphere completely inside the charge 
distribution. A reduced force is expected at the surface 
(see test particle at r2 in Fig.|4|. Both expectations are 
confirmed by the cumulative force, i.e., the integration 
over all charged shells, shown in Fig.[4]3. For this integra- 
tion we use Sr = (5^ — (r— r')^) / (2rr'), which guarantees 
that the "open" shells (see thick green/gray line in the 
sketch of Fig.[4| form a spherical correlation hole with 
radius ^ for \r — r^\ < 5. For all other shells it is (5t = 0. 
This Sr inserted into Eq. (|6| and integrated over all shells 
yields an explicit expression for the force in the presence 
of a correlation hole 

FR,6{r) = Qr/R^ for r < R-5 (7a) 

= An^sir) Qr/R^ for R-6 <r<R (7b) 

with the dimensionless attenuation factor 

{r^R-S)^ {2{r^R)S - 3{r-R)^^ 
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Indeed, the force from Eq. ^ increases linearly, charac- 
teristic for a homogeneously charged sphere and without 
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FIG. 5: (Color online) Energy spectra as obtained from an 
exploding cluster with 1000 particles for various cluster edges 
a (implemented with Fermi distribution function, see text), 
which is given in units of the cor rela tion hole radius 5. Energy 
is given in terms of E^^ cf. Eq. (3b). 



any effect of the correlation hole until the latter touches 
the surface from the inside (for the case shown in Fig.|4]3 
at r = 2R/i). When this happens the force grows more 
slowly than in the homogeneous case and even decreases 
still within the charged sphere reaching a maximal value 
at rmax = - ^/3 (for 8 < R). 

With the initial forces sculpturing the properties of 
the final ion-energy spectrum, we may even attempt to 
calculate this spectrum according to Eq. ([2| . We assume 
a self-similar expansion, i. e. scaling of all lengths in the 
system with the common factor 77. Consequently, the 
force ^ inherits the property Fr^ji^^jsiv^) = V~'^^R,6{^) 
from the Coulomb force. Thus, the final kinetic energy 
depends on the initial position similarly as in the case of 
a homogenous charge density through an integral along 
the similarity path with increment dr' = r dr] 
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In all (numerical and analytical) spectra presented we 
take into account a finite energy resolution 6E = 
For Eq. ([2| this means to replace S{x) with Kse{x) — 
exp( — (x/^£^)^)/y^(5£^. In Fig.jsjD one sees the spectrum 
for clusters with 1000 atoms in comparison to the one 
obtained with the mean-field approach including the cor- 
relation hole as just described. The deviations of both 



the numerical result (orange/gray-shaded area) and the 
correlation-hole result (solid line) with respect to the sim- 
ple mean-field result (dashed line) from Eq. ([3| at ener- 
gies E < E^ clearly reveals the importance of the corre- 
lation hole. 

However, given the analytically calculated force of 
Eq. ( 7b ) we also know that in order for the correla- 
tion hole to have an effect, the ion density must have 
a sharp edge, essentially fall to zero within the radius 
S of the correlation hole. Otherwise the forces of each 
shell in the presence of the correlation hole will com- 
pensate each other, see Fig. [4^. On the other hand, 
a soft edge gives rise to the mean-field peak through 
catching up of faster inner ions with the slower outer 
ones, as described in the introduction. We may quantita- 
tively describe the edge by a Fermi distribution Qa^nir) = 
2^a,i?/(l + exp((r— i?)/a)) with the softness parameter a 
and Qa^R the ion density at r=R which is determined 
through the integral Att j dr r'^ Qa^R(r)=N . In Fig.jsjthe 
energy spectrum of ions initially distributed according to 
a Fermi distribution Qa^Ri^) is shown for various values 
of a measured in terms of the correlation hole radius 5. 
One can see the cross over from the granularity to the 
mean-field peak near a = 5. 

We have seen that these two features are of very dif- 
ferent origin although both of them are surface effects 
and occur therefore on the rim of maximal energy in the 
spectra of Fig. [5] It should be possible in an experiment 
to reveal both phenomena and their cross over by sim- 
ply varying the pulse length of the ionizing light pulse 
as demonstrated in Fig.[l] A time-delayed probe pulse 
could reveal in addition the different temporal behavior 
of both peaks illustrated in Fig. [2] 

In the examples discussed here, the granularity of the 
system was quantified by the correlation hole which is the 
most common case for ground state matter. However, a 
little thought reveals, that ions (or electrons) randomly 
distributed over a finite volume, give rise to a granularity 
peak as well. It is much weaker and broader [T3| but still 
has the same reason: Forces on an ion from other ions 
inside and outside a virtual shell do not compensate each 
other near the edge of the sample where the radius of the 
shell is defined by the distance of the ion to the charge 
center of the sample. 

Acknowledgement — We gratefully acknowledge sup- 
port from CORINF, a Marie Curie ITN of the European 
Union, Grant Agreement No. 264951. 



[1] P. Fulde, Electron Correlations in Molecules and Solids 

(Springer, Heidelberg, 1995). 
[2] B. W. J. McNeil and N. R. Thompson, Nat. Photon. 4, 

814 (2010). 

[3] C. Gnodtke, U. Saalmann, and J.M. Rost, Phys. Rev. 

Lett. 108, 175003 (2012). 
[4] K. Nishihara, H. Amitani, M. Murakami, S. V. Bulanov, 



and T. Z. Esirkepov, Nucl. Inst. Meth. Phys. Res. A 464, 
98 (2001). 

[5] A. E. Kaplan, B. Y. Dubetsky, and P. L. Shkolnikov, 

Phys. Rev. Lett. 91, 143401 (2003). 
[6] Md. R. Islam, U. Saalmann, and J. M. Rost, Phys. Rev. 

A 73, 041201 (R) (2006). 
[7] C. Gnodtke, U. Saalmann, and J. M. Rost, New J. Phys. 



5 



13, 013028 (2011). 
[8] V. F. Kovalev and V. Y. Bychenkov, JETP 101, 212 
(2005), V. Y. Bychenkov and V. F. Kovalev, Plas. Phys. 
Rep. 31, 178 (2005). 
[9] F. Peano, R. A. Fonseca, and L. O. Silva, Phys. Rev. 
Lett. 94, 033401 (2005). 
[10] R. Garcia-Pelayo, J. Phys. A 38, 3475 (2005). 
[11] LJ clusters imply that the ground state configuration of 
those clusters is obtained by minimizing the LJ interac- 



tion between the atoms. 

[12] The ions have been placed radially on equal- volume shells 
guaranteeing that the resulting granular density approx- 
imates the homogeneous density as close as possible. The 
angular positions of the ions are relaxed to minimize their 
interaction energy given by mutual LJ potentials. 

[13] U. Saalmann, A. Mikaberidze and J. M. Rost, unpub- 
lished 



